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Abstract 

Vacuum expectation values (VEVs) of the current densities for charged scalar and Dirac 
spinor fields are investigated in (D + l)-diniensional de Sitter (dS) spacetime with toroidally 
compactified spatial dimensions. Along compact dimensions we impose quasiperiodicity con- 
ditions with arbitrary phases. In addition, the presence of a classical constant gauge field 
is assumed. The VEVs of the charge density and of the components for the current density 
along noncompact dimensions vanish. The gauge field leads to Aharonov-Bohm-like oscilla- 
tions of the components along compact dimensions as functions of the magnetic flux. For 
small values of the comoving length of a compact dimension, compared with the dS cur- 
vature scale, the current density is related to the corresponding current in the Minkowski 
spacetime by a conformal relation. For large values of the comoving length and for a scalar 
field, depending on the mass of the field, two different regimes are realized with the mono- 
tonic and oscillatory damping of the current density. For a massive spinor field, the decay 
of the current density is always oscillatory. In supersymmetric models on the background 
of Minkowski spacetime with equal number of scalar and fermionic degrees of freedom and 
with the same phases in the periodicity conditions, the total current density vanishes due to 
the cancellation between the scalar and fermionic parts. The background gravitational field 
modifies the current densities for scalar and fermionic fields in different ways and for massive 
fields there is no cancellation in the dS spacetime. 



PACS numbers: 04.62.-hv, 04.50.-h, ll.lO.Kk, 04.20. Gz 



1 Introduction 

The investigation of quantum effects in fixed gravitational backgrounds is among the most in- 
teresting topics in quantum field theory (for reviews see [I]). These effects may have important 
implications in black hole physics and in cosmology. The presence of the gravitational field, in 
general, reduces the number of symmetries and exact results can be obtained for highly sym- 
metric backgrounds only. A better understanding of physical effects in these backgrounds could 
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serve as a handle to deal with more complicated geometries. In this context, the investiga- 
tion of quantum field theoretical effects in de Sitter (dS) spacetime is of special interest. In 
inflationary scenario an approximately dS spacetime is employed to solve the problem of initial 
conditions in standard cosmology [2]. During an inflationary epoch in the early Universe the 
quantum fluctuations play a crucial role in the generation of cosmic structures from inflation. 
More recently observations of high-redshift Type la supernovae, galaxy clusters and cosmic mi- 
crowave background [3] indicate that at the present epoch the Universe is accelerating and is 
well approximated by a model with a positive cosmological constant as a dominant source for 
the expansion. If the Universe would accelerate indefinitely, the standard cosmology would lead 
to an asymptotic dS universe. Therefore, the investigation of physical effects in dS spacetime is 
important for understanding both the early Universe and its future. An interesting topic which 
has received increasing attention is related to string-theoretical models of dS spacetime and in- 
flation. A number of constructions of metastable dS vacua were discussed within the framework 
of string theory (see, for instance, [1]). 

In recent years much attention has been paid to the possibility for the universe to have 
non-trivial topology. In particular, a number of fundamental physical theories are formulated 
in spacetimes with compact extra dimensions. The idea of compactified dimensions has been 
extensively used in supergravity and superstring theories. The models of a compact universe 
with non-trivial topology may play an important role by providing proper initial conditions for 
inflation The quantum creation of the universe having toroidal spatial topology is discussed 
in [6] and in Refs. [7J within the framework of supergravity theories. As it was argued in Refs. 
[1], there is no reason to believe that the version of dS spacetime which may emerge from string 
theory, will necessarily be the most familiar version with symmetry group 0(1,4) and there 
are many different topological spaces which can accept the dS metric locally. There are several 
reasons to expect that in string theory the most natural topology for the universe is that of a 
flat compact three-manifold. 

The non-trivial topology of the background space can have important physical implications 
in quantum field theory. The periodicity conditions imposed on fields along compact dimensions 
give rise to the modification of the spectrum for normal modes and, related to this, the expecta- 
tion values of physical observables are changed. A well known effect of this type, demonstrating 
the connection between quantum phenomena and global properties of spacetime, is the topolog- 
ical Casimir effect [9l [TO] . The Casimir energy of bulk fields induces a non-trivial potential for 
the compactification radius providing a stabilization mechanism for moduli fields and thereby 
fixing the effective gauge couplings. The Casimir effect has also been considered as a possible 
origin for the dark energy in both Kaluza-Klein-type and braneworld models [llj . 

One-loop quantum effects for scalar and Dirac spinor fields induced by the toroidal com- 
pactification of spatial dimensions in dS spacetime have been recently investigated in Refs. 
[T2| [T3| [Ml [15] . In these papers the vacuum expectation values (VEVs) for the field squared and 
the energy-momentum tensor are considered for untwisted and twisted fields. These quantities 
are among the most important local characteristics of the quantum vacuum and are closely 
related with the structure of spacetime through the theory of gravitation. For charged fields 
another important bilinear characteristic is the VEV of the current density. In addition to 
describing the physical structure of the quantum field at a given point, the current acts as 
the source in the Maxwell equations and plays an important role in modeling a self-consistent 
dynamics involving the electromagnetic field. 

In the present paper, we consider the combined effects of topology and the gravitational field 
on the VEVs of the current density for charged scalar and fermionic fields in the background 
of dS spacetime with an arbitrary number of toroidally compactified spatial dimensions. Along 
compact dimensions we impose generic quasiperiodic boundary conditions with arbitrary phases 
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which, as special cases, include the periodicity conditions for untwisted and twisted fields. In 
addition, the presence of a classical constant gauge field will be assumed. As it will be seen, 
this leads to Aharonov-Bohm-like effects on the VEV of the current density. The paper is 
organized as follows. In the next section we describe the background geometry and evaluate 
the Hadamard function for a complex scalar field assuming that the field is prepared in the 
Bunch-Davies vacuum state. By using this function, in section [3l we investigate the VEV of 
the current density. The current density for the Dirac spinor field is discussed in section [H The 
main results are summarized in section [5l 

2 Hadamard function for a scalar field 

In the present paper the background geometry is described by the {D + l)-dimensional dS 
line-element in planar coordinates: 

ds^ = dt'-e^''^Y.^dx'f. (1) 
1=1 

The corresponding metric tensor is the maximally symmetric vacuum solution of Einstein equa- 
tions with the cosmological constant A = D{D — l)a~^/2. For the coordinates Xp = (x^, 
one has — cxd < < oo, I = 1, . . . ,p, and we assume that the coordinates Xg = [x'^^^, ...,x^), 
with q = D — p, are compactified to circles with the lengths Lf. < x^ < Li, I = p + 1, . . . , D. 
Hence, we consider the spatial topology x (S^)^. In the discussion below it will be convenient, 
in addition to the time coordinate t, to use the conformal time r = — ae~*/" with — oo < r < 0. 
For this coordinate the line element takes a conformally flat form with the conformal factor 
(a/r)^: ds'^ = [a/r)^ -q^ydx^^dx^ ^ where rj^y is the Minkowskian metric tensor. 

Firstly, we consider a complex scalar field ip{x) with a curvature coupling parameter ^, in 
the presence of a classical abelian gauge field A^. The corresponding field equation has the form 

{D^D'' + m' + iR) ^{x) = 0, (2) 

where = + ieA^ is the gauge-covariant derivative and e is the coupling between the 
scalar and gauge fields. For the scalar curvature in ([2]) one has R = D{D + 1)/q!^. In the 
most important special cases of minimally and conformally coupled scalars one has = and 
^ = (L> — 1)/(4L)), respectively. Since the background space is multiply-connected, in addition 
to the field equation one should specify the periodicity conditions along compact dimensions. 
We will assume generic quasiperiodic boundary condition 

(/?(t,Xp,Xg + Liei) = e'"'ip{t,Xp,Xg), (3) 

with constant phases a; and with being the unit vector along the dimension x', / = p+1, D 
(for physical effects of phases in periodicity conditions for fields in multiply-connected spaces 
see also [16]). The results below will be periodic functions of a/ with the period equal to 2tt. 
The special cases of untwisted and twisted scalar fields (periodic and antiperiodic boundary 
conditions), most frequently discussed in the literature, correspond to a; = and a/ = vr, 
respectively. 

Here we will consider the simplest configuration of the gauge field with = const. In this 
case, the gauge field can be excluded from the field equation by the gauge transformation 

A'^ = A^ + d^u, ip'ix) = e-''^^{x), CO = -A^x''. (4) 
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In the new gauge one has A'^ = and D'^ = V^. Now the quasiperiodicity condition takes the 
form 

(/?'(t,Xp,Xg + Liei) = e'"V'(*,Xp,Xg), (5) 

where 

ai = ai + eAiLi. (6) 

As it is seen, the phases in the periodicity conditions and the value of the gauge field are 
related to each other through a gauge transformation. In what follows, the physical results will 
depend on the phases in the periodicity conditions and on the gauge potential in the form of 
the combination dU]). Note that if A is the spatial vector with D components corresponding to 
the spacetime vector A^, then A; = —Ai. Although the corresponding field strength vanishes, 
a constant gauge field shifts the phases in the periodicity conditions along compact dimensions. 
As it will be seen below, this leads to the Aharonov-Bohm-like effects on the current density 
for charged fields. Note that in ([6]) the shift due to the gauge field may be written in the form 
eAiLi = —eAiLi = —2tt^i/^q, where $0 = 27r/e is the fiux quantum and is the flux enclosed 
by the circle corresponding to the lih. compact dimension. In what follows we will work in terms 
of the gauge transformed held (p'{x) omitting the prime. 

The physical quantity we are interested in is the current density for a charged scalar field. 
The corresponding operator is given by the standard expression: 

jf,{x) = ie[(p'^{x)D^(p{x) - {Df,(p'^{x))(p{x)]. (7) 

Note that in the gauge under consideration one has Df^ip = V^(/? = df^ip. The VEV of the current 
density, {jn{x)) = (0| |0), is expressed in terms of the Hadamard function 

G(i)(x, x') = (0| + ^"-{x'Mx) |0) , (8) 

where |0) stands for the vacuum state. The corresponding formula reads: 

0-^(x)) = ^elim(a^-a;)G«(x,x'). (9) 

By expanding the field operator in terms of a complete set {ip^a^\x)} of solutions to the field 
equation, obeying the periodicity conditions and using the standard commutation relations 
for the annihilation and creation operators, the following formula is obtained for the Hadamard 
function: 

G^'Hx,x') = Y.Y. V>i'\x)^^:>{x'). (10) 

Here includes the summation over the discrete components of the collective index a and the 
integration over the continuous ones (for the specification of the collective index in the geometry 
under consideration see below). 

Because of the plane symmetry of the problem, the dependence of the scalar mode functions 
on spatial coordinates can be taken in the standard exponential form, e*'' "^, with x = (Xp,Xg) 
and k = (kp,kg). From the field equation it can be seen that the general solution for the time- 
dependent part of the mode functions is a linear combination of the functions r]^/^Hi^\kri) and 
where /c = Ikl, 

r] = \T\= ae"*/", (11) 
(1 2) 

and the order of the Hankel functions ' (z) is related to the mass of the field by 

u = [1)2/4 - D{D + 1)C - m^Q^] . (12) 
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This parameter is either real and nonnegative or purely imaginary. 

An important step in formulating quantum field theory in a curved spacetime is the choice 
of vacuum. Different choices of the coefficients in the linear combination of the Hankel functions 
correspond to different choices of the vacuum state. dS spacetime has maximal symmetry and 
it is natural to choose a vacuum state having the same symmetry. It is well known that in dS 
spacetime there exists a one-parameter family of maximally symmetric vacuum states, which 
have been dubbed the a- vacua (see Ref. |17j and references therein). Among them the Bunch- 
Davies vacuum [18] is the only one with the Hadamard singularity structure. 

Here we assume that the scalar field is prepared in the Bunch-Davies vacuum. The mode 
functions realizing this state read: 

In these expressions, for the components of the momentum along uncompactified dimensions 
one has — oo < ki < +00, I = 1, . . . ,p, and the eigenvalues of the components along the compact 
dimensions are quantized by the periodicity conditions dS]): 

ki = {2TTni+ai)/Li, = 0, ±1, ±2, . . . , (14) 

for I = p + 1, D. Hence, the mode functions are specified by the set of quantum numbers 
a = (kp, rig), where = (up+i, no). 

The coefficients C^^'' in (jl3p are determined from the orthonormalization condition 



ip'^^\x)dt'f^P*{x) - ^^^?*{x)dt^i'\x)\ = i6^.'5ss', (15) 

with (5o-o-' = (^(kp — ^p) YliLp^i Snin[ ™d g being the determinant of the metric tensor corre- 
sponding to the line-element ([T|). By using the Wronskian relation for the Hankel functions, one 
finds 

where the star stands for the complex conjugate and Vg = Lpj^i.-.L^ is the volume of the 
compact subspace. 

Substituting the functions into the mode-sum formula (fTUj) , for the Hadamard function 
we get the representation 



V ' > 2P+27rP-iKa^-i ] ^ ^ 

AH^}\kr^)Iif)(kri)^Ii^;^){kri)B^}\kri)\ (17) 
where Axp= Xp — Xp, Axg= Xg — x^. In ()17p and in the formulas below we use the notation 

+00 +00 

n ni=— 00 ?!;=— 00 

for n = (ni, . . . ,n;). Note that for a scalar field in dS spacetime without compact dimensions 
{j) = -D), the corresponding two-point functions contain infrared divergences for Rei' > D/2 
and in this region the Bunch-Davies vacuum is not a physically realizable state. In particular. 



gikq-Ax, 
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this is the case for a minimally coupled massless scalar field. These divergences come from the 
singular behavior of the Hankel functions at the origin. For the topology under consideration, 
assuming that X^^p+i af ^ and \ai\ < ir, the momentum k has a nonzero minimum value 

fc min = \lYlb=p+i '^f / Lf- this case the two-point function (fT7|) contains no infrared divergences 
and the Bunch-Davies vacuum is a physically realizable state for all values of the parameter u. 

For the further evaluation of the Hadamard function, given by (fT7|l . we apply to the series 
over nr,p+l<r<D, the Abel-Plana-type summation formula \19\ [20] 



27r °° 



/'OO 

^ E 9{kr)f{\kr\)= dz[g{z)+g{-z)]f{z 

^ n,,=-oo 



+./ dz[n^z)-fi-^z)] Y: (19) 



A=±l 



with defined by (jl4p . In the special case g{z) = 1 and a,- = this formula is reduced to the 
standard Abel-Plana formula (for generalizations of the Abel-Plana formula see [21]). For the 
series in (jl7p one has g{z) = e*^^^"^ and the function f{z) is given by the expression in the square 
brackets. It can be seen that the contribution of the first integral in the right-hand side of (I19p 
to the Hadamard function coincides with the corresponding Hadamard function for the topology 
^P+i ^ (^S'^Y"^ with the lengths of the compact dimensions (-^^p+i, • • • , -L^-i, -L^+i, . . . , L/j). 

In the part corresponding to the contribution of the second integral in (jlOp we use the 
expansion — 1 = Yl^=i After some intermediate calculations, the Hadamard function is 

presented in the form 



^ '^ (2vr)(^+^)/V,a^-i ^0 

oo 

x[l^,{r^z)K,{i^'z) + K,{r^z)I,{r^'z)\ e 



z 



/(p-i)/2(\/^^T^y (Ax^' + nLrf + |AXp|2) 
[(Ax^ + nL^)V|Axp|2](p-i)/2 

(r) 



(20) 



where Iu{y) and K^{y) are the modified Bessel functions, n^_i = (np+i, . . . , n.f._i, n^+i, . . . , n^), 
k^'^i = {kp+i, kr-i, kr+i, ...,kD) and 

4 = \^tli\'= E (2vrn, + dz)Vi^?. (21) 
Here and in what follows we use the notation 



Uiz) = z^'K^iz). (22) 

The representation (j20p of the Hadamard function is valid for Rez^ < 1. The n = term in 
this formula is the Hadamard function for the topology RP'^^ x (S*^)*?^^ with the lengths of the 
compact dimensions (i^p+i, . . . , Lr-i-, Lr+i, ■ ■ ■ , Ld). The remaining part vanishes in the limit 
Lr — )• oo and it is induced by the compactification of the rth dimension. The expression pop 
can be used for the investigation of the VEVs for various physical observables. Here we consider 
the VEV of the current density. 
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3 Current density for a scalar field 



Having the Hadamard function, we can evaluate the expectation value for the current density 
by using the formula It can be seen that the VEVs of the charge density (/i = 0) and of the 
components of the current density along uncompactified dimensions vanish: 

{j^) =0, fi = 0,l,...,p. (23) 

The latter is a direct consequence of the problem symmetry. For a noncompact dimension 
the problem is symmetric under the reflection — )• —x^. The presence of the nonzero current 
density along would break this symmetry. 

For the VEV of the current density along the rth compact dimension, from ([9]) and (j20p we 
find 



X E E /(p-M)/2(nL.v/i^). (24) 

n=l (r) 
"q-1 



This VEV is an even periodic function of the phases ai, I ^ r, with the period 27r, and it is an 
odd periodic function of the phase otr with the same period. In particular, the current density is 
a periodic function of the magnetic flux with the period of the flux quantum. In the absence of 
the gauge field the VEV of the current density vanishes for special cases of twisted {otr = n) and 
untwisted (a^ = 0) fields. In the special case of a single compact dimension one has p = D — 1, 
g = 1, and the general formula takes the form 

= ^''"III+Cd 1 r ' [^-(^^) + ^'^(^^^l ^'^^^'^ E '^^^fD/2{nLrz). (25) 
(2tt) Lr JO 



From the covariant conservation equation (j'^) = it follows that the charge fiux through 
the {D — l)-dimensional spatial hypersurface x'" = const is determined by the quantity rir {j'^), 
where rir = \/\grr\ = Oi/rj is the normal to the hypersurface. Now, from (j24p we see that this 
quantity depends on the time coordinate and on the lengths of the compact dimensions in the 
form of the ratio Li/ij: 

rir if) = a-^/(Lp+i/7?, . . . , Ld/v)- (26) 

By taking into account that L^"^^ = aLi/rj is the proper length of the compact dimension, we 
see that the ratio Li/rj is the proper length of the Ith compact dimension measured in the units 
of the dS curvature scale a. 

For a conformally coupled massless scalar field one has m = 0, ^ = (D — 1) /(4L>) and, hence, 
v = 1/2. By taking into account the corresponding expressions for the modified Bessel function, 
we see that [I-^{x) + Iu{x)] Ky{x) = l/x. The integral in (f24j) is evaluated by using the formula 

duf^{a^v? + b^) = iy|/^+i/2(a6), (27) 
and for the VEV of the current density we get 

(27r) VqLr „=i ^ (,) 

"9-1 
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In this case, the current density is related to the corresponding result in the Minkowski spacetime 
(see Ref. [22]) with compact dimensions of the lengths (i^p+i, • • • ,Ld) by the equation (f) = 

in/a)^+^ (/)(^). 

Let us consider the Minkowskian limit for the general expression of the current density given 
by (j24p . This corresponds to the limit a — >• c« for a fixed value of t. In this limit one has 
I' ~ with /3 = ma ^ 1 and rj ^ a — t. By using the uniform asymptotic expansions for the 
modified Bessel functions for imaginary values of the order with large modulus (see, for example, 
Ref. [23]), we can see that for x < 1 one has 



[Iip{l3x) + I-ifiWx)] Ki^Wx) ~ ^ cos[2/3/(x)], (29) 



where 



i + \/r 



fix) = In y-—^ —j - Vl-x^. (30) 

In the case x > 1 the leading term is given by 

[I,pil3x) + I..pi/3x)] K,f,{l3x) ~ (31) 



From these expressions it follows that the dominant contribution to the current density in ([24 
comes from the integration range z > m. In this range, by using (j3ip and the integration formula 
([27|) . to the leading order we find 

4eLr^ sin(nar 



(/) - E ^ E /,/.«("i.v^rr^). (32) 

[ZTT) Vq 



(r) 



The expression in the right-hand side coincides with the VEV of the current density in the 
Minkowski bulk with toroidally compactified dimensions [22]. For this background geometry, 
the finite temperature corrections to the expectation value of the current density have been 
investigated in Ref. [22]. In this reference, we have also derived an alternative expression for 
the current density in the Minowskian bulk by using the zeta function approach. In the case of 
a massless field, the expression in the right-hand side of ([32]) reduces to the result which follows 
from (f28l) . 

Now we turn to the investigation of the current density in the asymptotic regions of the 
ratio Lr/rj. As we have mentioned before, this quantity is the ratio of the comoving length of 
the compact dimension to the curvature radius of dS spacetime. For small values of this ratio, 
Lr/i] ^ 1, we introduce in (124p a new integration variable y = L^z. By taking into account that 
for large values of x one has [I__j^(x) + Iu{x)] K^{x) w l/x, we find that to the leading order (j^) 
coincides with the corresponding result for a conformally coupled massless field (see (j28p ): 

if) « in/c^f^' O")^''^ , Lr/v « 1. (33) 

For a fixed value of the ratio Lr/a, the limit under consideration corresponds to the early stages 
of the cosmological expansion, t — )• — oo, and the current density behaves like exp[— (Z) + l)t/a]. 
Note that in the limit Lr/rj <C 1, the leading term in the VEV of the quantity rir (j^) coincides 
with the corresponding quantity in the Minkowskian bulk, if we replace the lengths of the 
compact dimensions L; by the proper lengths aLi/r]. In this limit, the effects induced by the 
curvature of the background spacetime are small. 

In the limit of large values of the proper length of the rth compact dimension, compared with 
the dS curvature scale, one has ry/Lr- <C 1. In this case, we introduce in ()24|) a new integration 
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variable y = L^z and then expand the integrand by using the formulas for the modified Bessel 
functions for small values of the argument. Two cases should be considered separately. For 
positive values of the parameter u, after the integration over y by using the formula from 
to the leading order we find 



if 



2^-(P-i)/2(ry/L,.) 



D-2u+2i 



7r(P+3)/2a«+iK 



n=l 



sm.{nar) 



E 

(r) 



f(j>+Z)/2-u{''^LrUr 



(34) 



n 
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for rj/Lr <^ 1. In the case of a conformally coupled massless scalar field = 1/2 and (fM|) 
reduces to the exact result given by ([28]) . For a fixed value of L^/a, the limit under consideration 
corresponds to late stages of the cosmological evolution, t — >• +oo, and the current density is 
suppressed by the factor exp[—{D — 21/ + 2)t/a\. Note that formula ([3^ also describes the 
asymptotic behavior for the current density in the strong curvature regime corresponding to 
small values of the parameter a. In the model with a single compact dimension {p = D — \) the 
asymptotic formula (j34l) takes the form 



2eT{u)T{D/2 + l-v)r] ^ sin(nQ 



7r^/2+ia^+i(L^/ry)S-2^+ 



lE 

n=l 



n 



D-2V+1 ' 



(35) 



The sum in the latter expression is expressed in terms of the Hurwitz zeta function ("(s,x). 

In the same limit, "q/L^ <C 1, and for pure imaginary values of the parameter i/, by taking 
into account that for small values of x one has 



[I-^,{x) + h{x)]K^{x) ^Re 



{2/xf^^\T{i\u\) 

r(i-^li^l) 



(36) 



for the current density we find the following asymptotic behavior 



(2vr)(P+3)/2^^^P- 
In this formula, Bg and (j)s are defined by the relation 



J cos[2|z^|t/a -|- 2\u\ ln(Lr/a) + 0^]. 



(37) 



2^Hr(i|z.|)£ 



n=l 



sin(nQr) 

j^+2-2i\v\ 



H hp 



-^)/2-i\v\{nLrUr 



(38) 



(r 
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Hence, in the case under consideration, at late stages of the cosmological evolution the current 
density is suppressed by the factor e"^^"*"^-**/" and the damping of the corresponding VEV has 
an oscillatory nature. 

For a fixed value of the time coordinate t, the condition ij/Lr <^ 1 corresponds to the length 
of the rth compact dimension much larger than the dS curvature scale, ^ a. Now for the 
values of the curvature scale a > one also has Lr 3> and the limit we have discussed 
corresponds to the length of the compact dimension much larger than the Compton wavelength 
of the field quanta. As it is seen from the formula (I35p . in this range the decay of the current 
density is a power-law as a function of mLr- This is in contrast to the case of the Minkowski bulk, 
where, as it can be seen from ()32p . the VEV of the current density is suppressed exponentially. 
For example, in the case of a single compact dimension {q = 1) one has 



if 



^(M) ^ 2em^sin(d^) 



(2vr) 



D/2 



m 



L,)^/2 
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for rriLr ^ 1. 

In figured] we display the quantity a^rir {j"^) /e as a function of the phase in the quasiperi- 
odicity condition for a conformally coupled scalar field in the D = 4 Kaluza-Klein-type model 
with a single compact dimension x"^ = (p = 3). The graphs are plotted for am = 0.25 and 
the numbers near the curves are the corresponding values for the ratio Lr/rj. Let us recall that 
the current density is a periodic function of ctr with the period equal 27r. For a field with the 
periodic boundary condition (a^ = 0) one has arjlix = — $r/*^o- 



2 
1 
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Q 
4S 
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Figure 1: The quantity a rir (j^) /e as a function of the phase in the quasiperiodicity condition 
for a conformally coupled scalar field in the D = 4 model with a single compact dimension. The 
graphs are plotted for am = 0.25 and the numbers near the curves are the corresponding values 
for the ratio Lr/rj. 

For the same model -D = 4, p = 3, in figure O we have plotted the quantity a^n^ (j^) for 
a conformally coupled scalar field versus the ratio Lr/rj. The latter is the proper length of the 
compact dimension measured in units of the dS curvature scale a. In numerical evaluations we 
have taken ar = it/2 and the numbers near the curves correspond to the values of ma. For small 
values of Lr/rj, from the asymptotic analysis given above it follows that {j'^) oc {Lr/r])~^ . 
As it has been explained before, depending on the parameter i^, for large values of Lr/rj two 
different regimes arise with monotonic (for positive values of z^) and oscillatory (for imaginary 
values of v) damping of the current density. Note that, for a conformally coupled field, is real 
for ma = 0.25 and imaginary for ma = 2, 3. In order to display the oscillatory behavior of the 
damping, on the right panel of figure [2] we have plotted the graph for ma = 3. The value of 
the ratio Lr/r] corresponding to the first zero of the current density decreases with increasing 
the value of ma. For the first two zeros of the current density one has Lr/ij = 3.88,8.29 and 
Lr/rj = 2.87,4.45 in the cases ma = 2 and ma = 3, respectively. 

4 Fermionic current 

In this section we consider the VEV of the current density for a fermionic field 'ip{x) on the 
background of {D + l)-dimensional dS spacetime with spatial topology x (S^)'^. As before, 
the line-element is taken in the form ([1]). Assuming the presence of a classical gauge field A^, 
the dynamics of the Dirac spinor field is governed by the equation 

ij^'Df,^ - = , = 5^ + + ieAf,. (39) 
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Figure 2: The current density for a conformally coupled scalar field in the model D = 4 and 
p = 3 as a function of the proper length of the compact dimension measured in units of the 
curvature scale a. The numbers near the curves are the values of the parameter ma, and for 
the phase in the periodicity condition we have taken = 7r/2. 

The Dirac matrices and the spin connection are expressed in terms of the flat-space Dirac 
matrices 7^"^ by the relations 

— ^(a) -p _ 1 (a) (b) (An.\ 
7 — ^(a) ' ' ^' ~ 4' ' ^{a)^{b)v;ii : I'iU; 

where the semicolon means the covariant derivative of vector fields and n are the tetrads 

(a) 

defined by ^^a)^'{b)^'^^ ~ 9^'^ ' ^^"^ background under consideration the tetrads can be taken 
in the form 

ef = 6% ef=e'l-b% a = 1,2,...,!). (41) 

In a (D + l)-dimensional spacetime the Dirac matrices are N x N matrices with = 2[(^"''^)/^] , 
where the square brackets denote the integer part of the enclosed expression. Here we take the 
fiat spacetime matrices in the Dirac representation 

7'°' - ( ; _°i ) • - ( -It ;° ) ■ (^^) 

with a = 1,2, . . . , D. 

We need also to specify the periodicity conditions for the fermionic field along compact 
dimensions. Here these conditions will be taken in the form 

V'(t,Xp,Xg + Lie,) = e'"'^/;(t,Xp,X5), (43) 

with constant phases a/. Although we have used the same notation for the phases in the 
boundary conditions for scalar and fermionic fields, of course, in general they can be different. 
The fermionic condensate and the VEV of the energy-momentum tensor in the special cases of 
untwisted (q; = 0) and twisted (a/ = vr) fermionic fields and in the absence of a gauge field have 
been investigated in [21 [15] . 

Assuming a constant gauge field A^, we can exclude it from the field equation by the gauge 
transformation A'^ = A^ -\- d^uj, 'ip'{x) = e~^^^'ip{x) with the function u defined in @. For the 
new field we have the quasiperiodicity condition 

ip'{t,Xp,Xq + Liei) = e^"'^'(t,xp,xg), (44) 
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with the new phase defined by (0). Similarly to the case for a scalar field, we will work in terms 
of the gauge transformed field omitting the prime. 

The VEV of the current density for the fermionic field can be expressed in terms of the 
two-point function sll\x,x') = (0|[V'r(aj), V's(2;')]|0)) where r and s are spinor indices. The 
expression for the VEV reads 

{fix)) ^ {0\fix)\0) = -^TiiYS^'\x,x)). (45) 

Expanding the field operator in terms of a complete set of solutions to the Dirac equation, 
{ip^\x)}, the following mode-sum formula can be obtained: 

O'"') = -| E E s^Pi'Hx)j>^iPi^\x). (46) 

In order to regularize the divergent expression in the right-hand side we will assume the presence 
of a cutoff function without writing it explicitly. The special form of this function will not be 
important for the further discussion. 

For the problem under consideration, the complete set of fermionic mode functions can be 
found in a way similar to that used in [T5] for special cases of twisted and untwisted fields. These 
mode function are presented in the form 

= 4+),(i^+i)/2,.k.x f ^l/L..^(^^)4^^ \ ^47) 

V Hl/2+iamikv)w'x J 

where n = k/A;, a = (fii, . . . , <t/)), and a = (k, x)- In (|47|l and (f48l) . Wx^\ x = 1, • • • , N/2, are 
one-column matrices having N/2 rows with the elements = and w\- = iw^j^\ For the 
normalization coefficients one has 

Similarly to the case of a scalar field, the eigenvalues for the components of the momentum along 
compact dimensions are given by the expressions (I14p . Note that for a massless field we have the 
standard conformal relation = [rj / a)^^'^'^^ ip^i^^^^^ between the mode functions defining 

the Bunch-Davies vacuum in dS spacetime and the corresponding mode functions ip^^'^^^^ in the 
Minkowski spacetime with spatial topology i?^ x (5'^)'^. 

Substituting the mode- functions into the mode-sum formula (j46p . it can be seen that the 
VEVs of the charge density and of the components along uncompactified dimensions vanish: 
{j^) = for fi = 0,1, ... ,p. For the component of the current density along the rth compact 
dimension we get the expression 



I- — 1/2— lam^ 1/2— jam^ — 1/2+jQm^ ''J' ^ ' 



where 



l/2+iam^ ''' 1/2— iam^ ''' — 1/2+jQm^ 

D 



kl+ (27rn, + dOV^'- (51) 

i=P+i 
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Further evaluation of the current density is similar to that for the case of the scalar field. 
Applying the summation formula ()19p to the series over rir in ()50p . one finds the expression 



27r)(P+3)/2a«+ 

oo 



— dzzRe{ [l-l/2-iamizri) + h/2+iam{zT])] (^??) } 

>/g Jo 



X E E fip-^minUV^^r), (52) 

n=l „W 
"9-1 

where is defined by (j2ip . Similarly to the case of the scalar field, the fermionic current density 
{j^) is an odd periodic function of the phase ctr with the period 27r and it is an even periodic 
function of the phases ai, I ^ r, with the same period. In the absence of the gauge field the 
current density vanishes for periodic and antiperiodic boundary conditions. Similarly to the 
scalar case, the quantity nr{j'^), which describes the charge flux through the hypersurface x^' = 
const, depends on r] and on the lengths of compact dimensions in the form of the combination 
Li/rj. The latter is the proper length of the compact dimension measured in units of a. 

For a massless fermionic field the modified Bessel functions are expressed in terms of the 
exponential function and, after the integration by using the formula (|27p . from (j52p we obtain 
the expression 

2e{r]/a)^+^N ^ sin(nd^) ^ 
^ = ~ (2vrW2+iKL? E E fp/2+iinLrUr). (53) 

^ ' y n=l (r) 

"9-1 

In this case one has the following relation between the fermionic and scalar current densities: 

(/)ferm = -N{f)sj2. (54) 

Here we assumed that the phases ai are the same for both fields. Note that in (j54p . {j^)sc is the 
current density for a complex scalar field which is equivalent to two real scalar fields. 

The transition to the Minkowskian limit is similar to that for the case of a scalar field. In 
this limit a — )• oo for a fixed t and r] ^ a — t. The dominant contribution to the integral in (|52p 
comes from the integration region z > m. In this region we have 

— 1/2— iam 

{zr]) + Il/2+iani{zr])] K 

1/2+ jam 

(zr?) ^ _ ^ . (55) 

ay z'^ — 

After the integration with the help of (j27p . to the leading order we get the current density in 
the Minkowski spacetime {f) ~ (j'')^^^ with 

oT«' = E ^ E /./..i(«x,V^Fw). (56) 

l^^J ^1 n=\ (r) 

This result coincides with the formula derived in |25j (in the notations of this reference a/ — )• 
— 27rd; and in [25] Ai = A; in the definition for ai ). The reason for the sign difference of a/ in 
the expression of d/ in ([6]) and in the corresponding expression of Ref. [25] is that, in the latter 
reference, for the evaluation of the VEV the negative-energy modes have been used, with the 
eigenvalues fc^"''^ = 27r(n/ + ai)/ Li instead of k\ = 2iT{ni — ai)/ Li. This means that, in fact, 
the formulas given in [25] are for the periodicity conditions ()43p with ai replaced by —2tt ai (see 
also the comment in [26]). Comparing (j56p with (j32p . we see that the relation (j54p between the 
fermionic and scalar current densities holds in Minkowski spacetime for massive fields as well. 
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Hence, in the super symmetric models on the background of the Minkowski spacetime with equal 
number of scalar and fermionic degrees of freedom the total current density vanishes. Note that 
this is not the case for the currents in the background of dS spacetime. 

Now let us consider the fermionic current density in the asymptotic regions of the ratio Ly/r]. 
For Lr/rj <^ 1 the dominant contribution to (j52p comes from the region with large values of ztj. 
By taking into account that for large x one has [l-i/2-iaTn{x) + Ii/2+iam.i^)] Ki/2+iaTn{^) ^ ^1^^ 
we see that to the leading order (j'^) coincides with the corresponding result for a massless 
fermionic field (j'') ~ {r\lot)^^^ (j^)^^^ ■ In this region we have the relation (j54p between the 
scalar and fermionic field currents. 

For small values of the ratio rj/Lr, the dominant contribution to the integral in ()52p comes 
from the range with small values of zr]. By making use of the expansions for the modified Bessel 
functions for small values of the arguments, to the leading order we find 

- 2P/2^(p+3)/2y^P + 2amln(L,/a) + 0/]. 



(57) 



Here, Bf and (j)f are defined by the relation 

oo 

i»"^r {1/2 + iam)Y, 



n=l 



sin(na,.) 



p/2-\-l—iam 



(58) 



Hence, unlike the scalar case, the damping of the fermionic current density for large values of the 
Lr/ij is always oscillatory. Another difference is that the oscillation frequency does not depend 
on the curvature scale of dS spacetime and is completely determined by the mass of the field 
quanta. 

Figure [3] presents the quantity a^rir (j^) /e for a fermionic field versus the phase in the 
quasiperiodicity condition for the D = 4 model with a single compact dimension (p = 3). The 
graphs are plotted for am = 0.25 and the numbers near the curves are the values for the ratio 
Lr/r]. 
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Figure 3: The current density as a function of the phase in the quasiperiodicity condition for a 
fermionic field in the D = A model with a single compact dimension. The graphs are plotted for 
am = 0.25 and the numbers near the curves are the corresponding values for the ratio L^/r]. 

The dependence of the current density on the proper length of the compact dimensions is 
displayed in figure H] for the same model = 4, p = 3. The graphs are plotted for the phase 
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in the periodicity condition corresponding to = 7r/2 and the numbers near the curves are 
the values of the parameter ma. For the fermionic field the damping of the current density for 
large values of Lr/rj is always oscillatory. In order to display this behavior, on the right panel of 
figure m we present the current density for the case ma = 3. As in the scalar case, the value of 
the ratio Lj./r] for the the first zero of the current density decreases with increasing ma. For the 
first two zeros we has Lr/r] = 3.04, 5.90 and Lr/ij = 2.51, 3.67 in the cases ma = 2 and ma = 3, 
respectively. 



0.00 




C5 




Figure 4: Fermionic current density as a function of the proper length of the compact dimension 
for separate values of the parameter ma (the numbers near the curves). The graphs are plotted 
for Or = it/2. 



5 Conclusion 

Among the most interesting consequences of the compactification of spatial dimensions is the 
appearance of nonzero expectation values for physical observables in the vacuum state of a 
quantum field. In the present paper we have considered the effect of the nontrivial topology 
on the VEVs of the current density for charged scalar and fermionic fields in the background 
of dS spacetime. In both cases we have assumed that the fields are prepared in the Bunch- 
Davies vacuum state. Along compact dimensions the quasiperiodicity conditions (j3]) and (I43p are 
imposed with general phases. In addition, we have assumed the presence of a constant classical 
abelian gauge filed. This leads to the Aharonov-Bohm-like effect on the current density. The 
physical quantities depend on the phases in the periodicity conditions and on the gauge potential 
in the form of the combination ([6]). The VEVs of the charge density and of the components of 
the current density along compact dimensions vanish. 

For the general case of toroidal spatial topology RP x {S^)'^, the VEVs for the scalar and 
fermionic current densities along the rth compact dimension are given by the expressions (I24p 
and ([52]) . respectively. These VEVs are even periodic functions of the phases a/, / / r, with 
the period 27r, and they are odd periodic functions of the phase with the same period. In 
particular, they exhibit Aharonov-Bohm-type oscillations as functions of the magnetic flux (with 
the period of the flux quantum) enclosed by the compact dimension. The current densities 
vanish in the special cases of untwisted and twisted flelds, when the gauge fleld is absent. 
We have explicitly checked the transition to the Minkowskian results, previously discussed in 
[22] [25| 126] . in the limit a — >• oo, for a fixed value of t. In this limit, if the phases in the 
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periodicity conditions and the masses for scalar and fermionic fields are the same, we have 
the simple relation (j54p between the scalar and fermioinic current densities. In particular, in 
supersymmetric models on the background of the Minkowski spacetime with equal number of 
scalar and fermionic degrees of freedom the total current density vanishes. The background 
gravitational field modifies the current densities for scalar and fermionic fields in different ways, 
and we have no similar cancellation in the dS spacetime. 

If the proper length of the compact dimension is much smaller than the dS curvature scale, 
the effects induced by gravity are small and the current densities are related to the corresponding 
results for massless fields in Minkowski spacetime with compact dimensions by the simple formula 
(j*") [r]/a)^^^ (j^)^^'^ ■ For a fixed value of the ratio L^/a, this limit corresponds to the early 
stages of the cosmological expansion (t — )• — oo), and the current densities behave as e"^^"*"^^*/". 
The effect of gravity on the VEVs is decisive for proper lengths of the compact dimensions larger 
than the dS curvature radius. In this limit, for the case of a scalar field, depending on the mass, 
two regimes are realized. For the first one, corresponding to real values of the parameter i/ 
defined by (jl2p . the leading term in the asymptotic expansion of the current density is given 
by (f34|) and the VEVs decay monotonically as {rj / Lr)^~'^'^~^'^ ■ For a fixed value of Lr/a, this 
corresponds to late stages of the cosmological evolution (t — )• +00), and the current density 
is suppressed by the factor ^-(^-"2^+2)1/ a ^ ^^iq second regime, realized for imaginary values 
of I', the behavior of the scalar current density is oscillatory, with the amplitude decaying as 
[t] / Lr)^'^'^ . The period of oscillations is given by 7ra/|z^|. For the fermionic field the oscillatory 
regime is realized only (see (j57p ) with the amplitude decaying as (rj/Lr)^^^. In the fermionic 
case the period of oscillations does not depend on the curvature radius and is given by 27r/m. 
Note that the decay rates for the amplitude in the cases of the scalar and fermionic currents 
are different. At late stages of the cosmological expansion, the current densities are suppressed 
by the factors e"'-^"''^^*/" and e^*-^^"*^^*/" for scalar and fermionic fields, respectively. For the 
values of the curvature scale a > m^^, if the length of the compact dimension is much larger 
than the Compton wavelength of the field quanta, mLr ^» 1, the decay of the current densities 
as functions of mLr is a power-law, for both scalar and fermionic fields. This is in contrast 
to the case of the Minkowski bulk, where the VEVs of the current densities are suppressed 
exponentially. 

The effects we have discussed can be applied to two types of models. For the first one D = 3 
and the results given here describe how the properties of the universe are changed by one-loop 
quantum effects induced by the compactness of spatial dimensions. Though the observational 
data constrain the size of possible compact dimensions to be larger than the horizon scale, at 
early stages of the cosmological expansion the lengths of compact dimensions were small and 
the effects considered here can be important. The second class of models, with p = 3 and 
D > 3, correspond to the universe with Kaluza-Klein-type extra dimensions. In these models, 
the current density along compact dimensions could be a source for cosmological magnetic fields. 
The gravitational analog of this electromagnetic effect was discussed in [llj, where it was shown 
that the topological Casimir energy can be considered as a possible origin for the dark energy. 
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